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^ ■ Abstract. 

ON , The chemical abundances in damped Lyman-a systems (DLAs) show 

' more than 2 orders of magnitude variation at a given epoch, possibly 

^ . because DLAs arise in a wide variety of host galaxies. This could signif- 

Q\ ' icantly bias estimates of chemical evolution. We explore the possibility 

<^> . that DLAs in which H2 absorption is detected may trace cosmological 

chemical evolution more reliably since they may comprise a narrower set 
of physical conditions. The 9 known H2 absorption systems support this 

I hypothesis: metallicity exhibits a faster, more well-defined evolution with 

■ redshift than in the general DLA population. The dust-depletion factor 

\ and, particularly, H2 molecular fraction also show rapid increases with 

O { decreasing redshift. We comment on possible observational selection ef- 

' \ fects which may bias this evolution. Larger samples of H2-bearing DLAs 

^ ■ are clearly required and may constrain evolution of the UV background 

' and DLA galaxy host type with redshift. 

a- 

> ' 

j_j ■ 1. Introduction 

. p.: 

Apart from providing independent supporting evidence of the big bang, the de- 
tection and subsequent study of cosmological chemical evolution provides the 
empirical details of galaxy formation and evolution. How primordial and pro- 
cessed gas is consumed by star formation, the dominant feedback processes and 
merging scenarios, may all contribute to the overall evolution of chemical abun- 
dances. One high-precision probe of this evolution is the spectroscopic study 
of damped Lyman-a systems (DLAs): absorbers with neutral hydrogen column 
densities A r (H 1) > 2 x 10 20 cm . Although these observations demonstrate that 
DLAs arise along lines of sight through distant galaxies, they do not directly dis- 
close details such as the galaxy's morphology, luminosity, mass or age. 

There exists substantial evidence that DLAs arise in a variety of galaxy 
types. At low-z (z a b s ^ 1-5), kinematic and Hi 21-cm absorption studies 
(e.g. Briggs et al. 1985; Briggs, de Bruyn, & Vermeulen 2001) suggest a signifi- 
cant contribution from spiral galaxies, a view supported at high-z by kinematic 
modelling and abundance studies (Prochaska & Wolfe 1997; Wolfe 2001). How- 
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ever, direct imaging at low-z (e.g. Le Brun et al. 1997) reveals that DLA hosts 
are a mix of irregulars, spirals and low surface-brightness galaxies (LSBs). A 
recent z = 0Hi 21-cm emission study (Ryan- Weber, Webster, k, Staveley-Smith 
2003) supports this. Boissier (2003) argue that the number of DLAs per red- 
shift interval and the iV(Hi) distribution imply that DLAs at z < 2 are a mix of 
spirals and LSBs whereas, at higher z, they are more likely to be dwarfs. This 
is supported by fitting of chemical evolution models to DLA metal abundances 
(e.g. Baker et al. 2000) and by recent 21-cm absorption searches at high redshift 
(Kanekar & Chengalur 2003). Though further work is clearly needed, the direct 
and indirect evidence for a 'mixed bag' of DLA hosts is already compelling. 

2. Selecting a less biased tracer of chemical evolution 

Evidence for an increase in DLA metallicities, [M/H], with cosmic time has 
emerged only gradually (Pettini et al. 1995; Lu et al. 1996; Vladilo et al. 2000; 
Kulkarni & Fall 2002; Prochaska et al. 2003), the latter reference using over 100 
DLAs to provide the strongest statistical evidence so far. Such large samples are 
required due to the huge scatter (~2dex) in [M/H] at a given epoch (see Fig. 1), 
a diversity expected given the variety of DLA hosts discussed above. However, 
the diversity in [M/H] could also significantly bias any estimate of chemical 
evolution, as could several observational selection effects (e.g. Hou et al. 2001). 

In Curran et al. (2003) we suggest that by selecting those DLAs in which H2 
absorption is detected, one may reduce or possibly avoid some of the biases be- 
setting the general DLA population. There are currently 9 confirmed H2-bearing 
DLAs (see below) and, typically, H2 is detected in only a few velocity compo- 
nents. These H2-bearing components seem distinct from the others, showing 
lower temperatures and higher dust depletion factors [M/Fe] (Petitjean et al. 
2002). Therefore, H2-bearing DLAs might be a less biased tracer of chemical 
evolution than the general DLA population since they may allow one to focus 
on a narrower range of physical conditions throughout cosmic time. 

2.1. Known H2-bearing DLAs 

H2 is detected in DLAs via the Lyman and Werner-band UV absorption lines 
which generally lie in the Hi Lyman-a forest. A compilation of results from H2 
searches in DLAs is given in table 8 of Ledoux, Petitjean, & Srianand (2003, 
hereafter L03), for which 7 DLAs have confirmed H2 detections and metallicity 
measurements: 0013-004 (z abs = 1.973), 0347-3819, 0405-443 (z abs = 2.595), 
0528-2505, 0551-366, 1232+082 and 1444+014. The DLA towards 0013-004 
(Petitjean, Srianand, & Ledoux 2002) comprises several absorbing components, 
of which 2 are dominant (their components c and d): iV(Hi) is measured at 
the mean redshift and we use the mean iV(H2) with error given by the range 
of iV(H 2 ) in L03. For the DLA towards 1232+082 we use the N(R 2 ) value and 
error from Srianand, Petitjean, & Ledoux (2000). 

We include two further DLAs: (i) The recent detection towards 0515—4414 
(Reimers et al. 2003) and (ii) 0000—2620, regarded as only a tentative detection 
by L03. However, the H2 identification has been carefully scrutinised by Lev- 
shakov et al. (2000, 2001) and relies on two H2 absorption features, the L(4-0)R1 
and W(2-0)Q(1) lines, the former appearing relatively clean from Lyman-a for- 
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Figure 1. Metallicity [M/H], dust-depletion [M/Fe], and molecular 
fraction (/) evolution for confirmed (filled circles) and tentative (hollow 
circle) H2 detections compared with values from the Prochaska et al. 
(2003) DLAs (small diamonds). The values adopted are discussed in 
the text. Linear fits to the H2-bearing DLAs (long dashes) and general 
DLAs (short dashes) were obtained by quadrature addition of a con- 
stant to the individual error bars such that = 1. Due to its much 
lower z a b s and potentially high photo-dissociation rate (Reimers et al. 
2003), we do not include 0515—4414 in these fits (but see Table 1). 



est blending. A tentative detection of H2 in 0841+1256 has also been reported 
(Petitjean et al. 2000), though confirmation requires future data and analyses. 

2.2. Results 

In Fig. 1 we plot against z a b s the metallicity [M/H], dust-depletion factor [M/Fe] 
and molecular fraction / = 2iV(H 2 )/[iV(Hi) + 2iV(H 2 )] for the 9 H 2 -bearing 
DLAs, comparing the former two quantities with DLAs in Prochaska et al. (2003) 
with Zn, Si, S or O metallicity over the relevant redshift range. [M/H] for the 
H 2 -bearing DLAs is [Zn/H], except for 0347-3819 and 1232+082 where it is 
[S/H] and [Si/H]. For 0515-4414, [Zn/Fe] is from de la Varga et al. (2000). 

The main results from Fig. 1 are summarized by the statistics in Table 
1: (i) [M/H], [M/Fe] and / for the H 2 -bearing DLAs are anti-correlated with 
£abs & t the 95% confidence level (i.e. more significant than for the general DLA 
population), (ii) [M/H] shows a steeper evolution with z a b s and a smaller scatter 
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Table 1. Statistics for correlations in Fig. 1. -P(r) is the probability 
of a chance correlation using a Kendall's-r test, B is the slope (with 
1 <r error) of the best fit, F is the F-test statistic (ratio of variances of 
two distributions) derived after subtracting the best fits from the data, 
and P(F) is the probability that F could be exceeded by chance alone. 
Since 0515—4414 is at much lower z a b s and may have a much higher 
photo-dissociation rate than the other H2-bearing DLAs (Reimers et al. 
2003), we present separate statistics for samples excluding and includ- 
ing 0515—4414. Note that this makes little difference. 
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about the slope than the general DLA population, (hi) [M/Fe] in ELi-bearing 
DLAs shows strong evolution with z a bs while the general DLAs show no evidence 
for evolution, and (iv) / ranges over 6 dex and shows a very steep evolution with 
z a bs- The new results (i)-(iii) support our hypothesis that ILi-bearing DLAs form 
a chemically distinct sub-class and may trace chemical evolution more reliably. 
The /-z a bs correlation was studied by L03 and is discussed further below. 

The [M/H]s and /s in Fig. 1 are measured using the total 7V(Hl) across 
the DLA profile and are not specific to the H2-bearing components. [M/Fe] is 
generally found to be uniform across most DLA profiles (Prochaska 2003), in- 
dicating that [M/H] should be uniform. However, the H2-bearing components 
typically have much higher [M/Fe] values than other components in the same 
DLA (e.g. L03 and 0347-383's [Si/Fe] profile in Prochaska 2003). These compo- 
nents usually dominate the non-refractory metal-line profiles and so, although 
[M/H] and / will be systematically underestimated, the effect will not be large. 
The fitted slopes in Fig. 1 are likely to be reasonably robust against this effect, 
but a larger sample and more detailed study is clearly required. 

What observational selection effects and biases could contribute to the steep 
/-z a bs evolution? Firstly, the sample is inhomogeneous since the spectra do 
not all have similar S/N and since the H2 detection methods and criteria were 
not uniform. Indeed, the weak H2 lines detected towards 0000—2620 are at the 
typical non-detection level (see fig. 16 in L03). The /-z a bs correlation is therefore 
only tentative. Secondly, the H2 detection limit will alter with z: equivalent 
widths increase but Lyman-a forest blending worsens with increasing z. Though 
this is an unlikely culprit for the 6 dex evolution observed, precise quantification 
of these competing effects requires numerical simulations. DLAs containing large 
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amounts of dust could suppress detection of their background quasars and may 
therefore be 'missing' from our sample (Fall, Pei, & McMahon 1989). However, 
since [M/H], [M/Fe] and / are positively correlated with each other (e.g. L03), 
this effect is likely to suppress, rather than create, the correlations in Fig. 1. A 
recent survey for DLAs towards radio-selected quasars (Ellison et al. 2001) also 
indicates that the number of such 'missing' quasars is likely to be small. 

3. Interpretations and conclusions 

Selecting those DLAs which exhibit H2 absorption may focus on systems with 
a narrower range of physical conditions than the DLA population as a whole. 
Tentative support for this conjecture lies in the steeper, tighter [M/H]-z a b s anti- 
correlation observed for the H2 systems studied here. Hou, Boissier, & Prantzos 
(2001) recently presented detailed chemical evolution models which give a slope 
for the [M/H]-2; a bs relation of ~ — 0.6dex. They correct this result for various 
observational biases to match the shallower slope observed for the general DLA 
population. However, the steep [M/H]-z a b s slope observed for H2-bearing sys- 
tems could be less affected by these biases and may avoid those introduced by 
sampling many different ISM gas phases. It might therefore be more compara- 
ble to the uncorrected slopes in the models. H2-selected DLAs are therefore a 
candidate for a less biased probe of chemical evolution. 

That there exists such a large range (~ 6 dex) in the values of / in Fig. 1 
may not be surprising: Schaye (2001) describes a [Zn/H] = —1 photo-ionization 
model for clouds in local hydrostatic equilibrium. For a representative inci- 
dent UV background flux and dust-to-metals ratio, the molecular fraction in 
this model shows a sudden increase of ~ 4 dex for only a small increase in the 
total hydrogen density. Therefore, the very steep /— £ a bs correlation could be 
achieved with a modest increase in dust content at lower redshifts, consistent 
with the observed [M/Fe]-z a b s anti-correlation. L03 also discuss how / might 
be very sensitive to local physical conditions. For example, within the Schaye 
model, one expects an anti-correlation between / and the intensity of the UV 
background. However, the behaviour of the UV background flux with redshift 
over the range z = 1-3 is still a matter of considerable uncertainty. The strong 
decrease in / at high redshift may also be consistent with recent H 1 21-cm 
absorption measurements in DLAs (Kanekar & Chengalur 2003), where a gener- 
ally higher spin/excitation temperature is found at z > 2.5. With an increased 
sample size and more detailed analyses, the /-£ a bs anti-correlation, if real, may 
provide complementary constraints on these problems. 
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